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The propagation of a spherical wave through a medium containing anisotropic random 
irregularities is considered. The formulation follows closely that of Karavainikov. The 
mean square deviations and the correlation functions of the phase and of the logarithmic 
amplitude are derived by assuming that the autocorrelation function of the dielectric con- 
stant is Gaussian with ellipsoidal symmetry. This form of autocorrelation function is 
chosen because experimentally it has been found that irregularities in the ionosphere at 
heights of 300 kilometers are cigar-shaped with approximately Gaussian ellipsoidal sym- 
metry. Methods of generahzing to other correlation functions are also indicated. 

The theory is applied to a problem which is of particular interest in the investigation 
of ionospheric irregularities by means of radio transmissions from satellites. Specifically 
the d(ipendence of scintillation (phase and logarithmic ampUtude) on the height of thi; trans- 
mitter above a slab of irregularities and the dependence of the autocorrelation functions are , 
investigated. The theory explains that the scintillation index should be relatively insensi- 
tive to the zenith angle of the satellite position for a temperate latitude station, in agreement 
with the present preliminary observations. The theory also suggests that the observation 
of phase scintillation may yield information concerning the physical location and the thick- 
ness of the region of irregularities. 

1. Introduction 

The study of wave propagation in a randoJii Jiiediuni is not only tlieoretically interesting, 
but also has wide applications in many brandies of physics, including acoustics, ionospheric 
and atmospheric physics, and astronomy. The random irregularities would scatter waves 
and hence cause wave interference. As a residt the amplitude and the phase of the wave 
fluctuate. The study of these fluctuations can yield certain statistical properties of the medium 
although, unfortunately, it usually cannot be done uniquely. 

The early foundation of the theoretical work was laid by Bergman [1946], who used ray 
theory, and Pekeris [1947], who studied the problem of scattering. Since then, many papers 
have appeared in the literature and Keller [1960] has even classified the mathematical methods 
used into '^honest" and '^dishonest'' methods. More recently, Silverman has translated the 
books of Chernov [1961] and Tatarski [1961], which largely represent the Russian contributions 
to the field. The present study is motivated by the observation of scintillation of radio signals 
received on the ground from an orbiting satellite-transmitter [Yeh and Swenson, 1959]. It 
is hoped that this paper will demonstrate that new information can be obtained by studying 
these scintillations and their relation to ionospheric irregularities. 

The analysis in this report is carried out with the application to wave propagation through 
ionospheric irregularities in mind. The region of irregularities is found to be above about 300 
km [Swenson and Yeh, 1961]. There is indication that the region may even be as high as 
1,000 km in the auroral zone [Easier and DeWitt, 1962]. Its occurrence as evidenced by the 
radio-star scintillation studies is correlated with the spread F phenomenon observed in iono- 
grams [Lawrence, Jespersen, and Lamb, 1961]. These irregularities mainly appear near mid- 
night [Booker, 1958; Yeh and Swenson, 1959]. It should be pointed out that tropospheric 
irregularities have negligible effects in the present problem. (For an excellent review article 
on tropospheric scattering, see Wheelon, 1959.) 

In this report it is assumed that the inductive permeability of the medium is that of free 
space, but the dielectric constant is a random variable of the position. The medium is char- 
acterized by the refractive index, 
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N{x)=<7,y[l+e^^ix)], (1> 

where <C^> is the smoothed average background value and is assumed to be a constant, e is a 
small constant parameter and /i is a random variable of position. Accordingly, the average 
value of jLt is zero. It shall be assumed that the correlation function of jjl is only a function of 
the difference of coordinates. In a medium of weak random irregularities the percentage change 
of the refractive index from its average value is small. This means that e is small and may be 
used as an expansion parameter. 

The problem of interest is essentially that of the diffraction of a spherical wave by a slab 
of intervening medium consisting of random irregularities. These irregularities are assumed 
to be so weak that only first-order scattering is of importance (Born approximation). In such 
a problem many characteristic lengths are involved. A brief discussion on these lengths is 
carried out in section 2. The formulation of the problem in section 3 follows closeW that of 
Karavainikov [1951]. We start from the scalar Helmholtz wave equation and apply the method 
of perturbation. The solution is expressible in terms of a series in which the zeroth order term 
is the solution in the absence of the irregularities (i.e., the familiar spherical wave solution). 
The first-order term shows the effect of the wave being scattered once, and the remaining 
higher order term shows the effect of multiple scattering. In a medium whose property varies 
sufficiently slowly over one wavelength these higher order Born solutions may be neglected. 

For purposes of obtaining explicit expressions the autocorrelation function of the fluctua- 
tions of the dielectric constant must be assumed. In the ionosphere the exact form of this 
correlation function has not yet been accurately determined. Experiments have suggested 
that the irregularities have a prolate spheroidal shape [Spencer, 1955]. The ratio of the major 
to minor axis is approximately 5. The length of the minor axis is of the order of 1 km. Hence 
one of the most important points that should be emphasized here is that these irregularities 
are highly anisotropic. For convenience the autocorrelation function is assumed to be 

p/x)=exp-(xyil+yyt,+z'lll), (2) 

in sections 4, 5, and 6, in which expressions for the mean square values, and the transverse and 
longitudinal correlations are obtained, respectively. When the correlation function departs 
from (2) some of these results are still useful. A method of generalization is presented in 
section 7. 

Applications of the theory to ionospheric studies are considered in section 8. For the 
numerical values chosen we find the scintillation of the logarithmic amplitude increases mono- 
tonically as a function of the height of the transmitter above the region of irregularities while 
that of the phase increases sharply at first when the transmitter is in the region and reaches a 
maximum near the top of the region and decays slightly above the region. Sample representa- 
tive curves for the autocorrelation functions are also plotted. The existing experimental 
results are still very sparse, but rough agreement between the theory and the experiments has 
been found. 

2. Characteristic Lengths 

The geometry of the problem is shown in figure 1 where a, b, and c respectively represent 
the distance between the transmitter and the top of the slab, the thickness of the slab, and the 
distance between the receiver and the bottom of the slab. The 2- axis is in the direction AB. 
In case the transmitter is inside the slab we shall ignore the backscattering and hence the 
effective thickness is reduced accordingly. As a approaches infinity the problem reduces to 
that of scattering from a uniform plane wave. 

It is convenient to normalize the lengths by the wave number k defined by, 

k==<^<n>/c, (3) 

whsre CO is the angular frequency of the wave and c the velocity of light in free space. Unless 
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otherwise specified it shall bo understood that all lengths are measured in this dunensionless 
unit. 

In the analysis of tlie scatteiiiiii- problem shown in figure 1 a few geometric parameters are 
of importance. TluM-c^fore, a bi-i(^r discussion of them is in order. 

A (TRANSMITTER) 



n 



Oo<^o%0 



SLAB OF RANDOM^ 



IRREGULARITIES, (^0 , 



o 








B (RECEIVER) 
Figure 1. Geometry of the problem. 



Z-y' 




RECEIVER 




RECEIVER 2 



Figure 2, Scattering from a single irregularitij. 



Figure 3. Magnification parameter. 



Consider tiie case of scattering froiu a single irregularity of size I as illusti'ated in iigure 2. 
When tlie condition, 

r,+r,-z«ll2, (4) 

is satisfied the phase change of the scattered waves at the receiver B introduced by the path 
difference is small. The scattered waves form almost a parallel beam and the fluctuation of the 
phase is predominately due to the change in refractive index along the irregularity. This case 
is usually referred to as Fresnel diffraction. On the other hand, when the converse of (4) is 
true it is then in the regime of Fraunhofer diffraction. 
For convenience define a wave parameter, 



In terms of this parameter the following equivalent statements can be made: 

D^y>l region of Fresnel diffraction, and 
Z^<C<Cl region of Fraunhofer diffraction. 



(5) 



(6) 



For anisotropic irregularities, since the scale of the irregularity depends on the direction it is 
then necessary to define several wave parameters. These are discussed in a later section. 

Another parameter that enters into tlie study of correlation functions can be seen as 
follows. In optics an object of size / has its shadow given by Xo=2l/y^ (see fig. 3). For this 
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reason one can define a ^^magnification^' parameter, 

M=xlxQ=xy'lzl, (7) 

Therefore, in the language of geometric optics, for a given 7' we would expect that the correlation 
function of the diffraction pattern on the ground has appreciable value only for values of x such 
that M is of the order of unity. 

3- Formulation 

In a slowly varying medium in which the property of the medium changes little in one 
wavelength the Helmholtz wave equation is approximately satisfied [Stratton, 1941], In the 
normalized unit, it is written as (time dependence exp {joit) has been assumed), 

VV+[l+6M(a:)]V=0. (8) 

When the medium is homogeneous (i.e., €=0) a particular solution of (8) for an isotropic point 
source yields the well-known spherical wave function, 

yPo= (kAo/r) exp (-jr) . (9) 

The problem is formulated by assuming the solution of the form, 

Hx)=Mx) exp -jlect>i{x) + e'<t>2ix)+ . . . ]. (10) 

Substitute (10) into the wave equation and equate coeflScients of equal powers in e. The 
result is a chain of equations. With suitable substitutions [Karavainikov, 1957] each of these 
equations can be transformed into an inhomogeneous wave equation and solved in the regular 
manner. When transformed back, the <^'s in (10) are given by 

<^i= W^^)^ ^ exp [-Jir'+R-r)wt', (11) 

<t>2=iJrMl f'^'^'y-J^^'t'iix')? ^^p [_j(^,^+E-r)]d^x', (12) 

where R is the distance connecting the source point (in this case the scatterer) and the observing 
field point (see fig. 1), i.e., 

R=\x-x'\, (13) 

The expressions given by (11) and (12) are physically reasonable. Each element of the 

transmitted spherical wave exp (—jr')lr' is scattered at x^ and reaches the receiver as another 
spherical wave exp (—jR)/R. Proper weighting functions have been introduced to take 
care of the amount of scattering. In (11) the weighting function is proportional to the 
percentage change of the refractive index. The weighting function in (12) has two terms; both 
are of higher order. Because of the appearance of the gradient term in (12), the expression 
for 02 is rather complex. From here on it is assumed that the irregularities are so weak that 
4)2 can be neglected. Since, as shown in (14), 01 is related to the phase and logarithmic ampli- 
tude, the neglect of the gradient term in (12) is equivalent to the assumption that the change 
of these quantities be small in one wavelength. 

It is convenient to rewrite the wave function (10) as, 

^l^(x) = (kAo/r) exp l-j{r+e<l>i)]=[kA(x)/r] exp ^j[r+Q(x)]. (14) 
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Note that when written this w^ay A represents the ampHtude of the wave and Q the phase 
departure. Substitute (11) into (14) and expand the exponential to the second order. For 
cases of scattering- from a locahzed region the resulting expression can be shown to be identical 
to that derived by Booker and Gordon [1950], except for the factor that takes into account 
the polarization of the incident wave with respect to the direction of scattering. Since the 
present analysis is primarily concerned with the scattering in the forward direction this factor 
reduces to 1. 

Put (11) into (14) and equate separately the real and the imaginary parts, 

Q{t) = {r^/2T)j^ ^ sin (r' +R-r)d4' , (15) 

log (£)=irel2^)j^ ^ cos (r'+R-r)(fx' . (16) 

It is known that a wave cannot be used to reveal structures much less than one wave- 
length [Ratcliffe, 1956]. The case of general interest is that in which the characteristic scale 
of the irregularities is much larger than one wavelength. In the case of the ionosphere the phys- 
ical size of I is roughly 1 km which corresponds to / (in the normalized unit) of the order of 400 
at 20 Mc/s. In our normalized unit this is equivalent to the condition, 

|Z»1. (17) 

Because of this condition the contribution to the integrals (15) and (16) comes predominately 
from the scatterers in the neighborhood of the straight line coimecting the transmitter and the 
receiver. Therefore, certain approximations may be used to facilitate the evaluation of these 
integrals. 

Define the normalized autocorrelation function by, 

p.(^)=<m(xi)m(?2)>/<m'>. (18) 

Here the symbol O is used to denote the spatial average and the x is the relative coordinates 
defined b^', 

a:=a:2— ici. (19) 

It shall be assumed that the ergodic hypothesis is valid so that the ensemble average is 
equivalent to the spatial average. It is seen immediately from (15) and (16) that the mean 
value of the phase departure and the logarithmic amplitude vanish, 

<Q>=0, (20) 

<S>=0, (21) 

where for simplicity S is written for log (A/Aq). 

The correlation functions can be represented formally by, 

^/i/^\/i/~^w /-^ 2^ 2^iA 2\ C C sin (ri+i?i— n) sin (r2+S2— ra) /""/n 73"1 «"^. 

(22) 
and 

^o/"^^o/"'^\ / 2^ 2^iA 2^ r r cos (r[+Ri—ri) cos (r2+jB 2— ^2) /^,n 73"", 73"^, 
<S(xi)S(x2)>={rir2e^<fi'>I^Tr') -Tp -7^ p^{x')d^x[(Px2. 

Jcj J V2 '1^1 "2-^2 

(23) 
625 



When X2=a:i, the two receivers coincide and the above expressions reduce to mean square 
values. 

4. Mean Square Values 

-^ -> 
For the purpose of computing the mean square values let 2:2= ^^1= (0,0, 3) in (22) and (23). 
It is convenient to make the following approximations. When the quantity appears in the 
phase let (see fig. 1), 

and when the quantities appear in the denominator 

r'^z' ,E^z-z'. (25) 

Define, 

h-«Q'> + <S'»le'<y?>, (26) 

and 

/2=7r« Q'> -<S'»le'<y?> . (27) 

The expressions for these integrals can be found b}^ the use of (22) and (23), 



where the new symbols are defined by 

^[=z[{z-2[)/z, and ^2=^2(^ — ^2)/^. (30 

In order to integrate (28) and (29) it is desirable to make the coordinate transfonnation 
by introducing the relative coordinates, 

X'-=^X2 — X[, y'=y'^ — y[^ ^' = 2^ — sj, (31) 

and the center of mass coordinates 

a' = (x[ + x:,)/2, p^ = {y[+y',)/2, Y = {z[ + z',)l2. (32) 

Integrate (28) and (29) with respect to a^ and I3\ resulting in, 

i,= r^' r r r ^-^^^ sin ^^z^^ dx^dy^dz^d^, m 

(* (l-\-l) f* I) /'CD /^ CO /'\ ' \ f 

'-J,.., J...-. I.-. L-. 2Wffl "" 2T?iTk) ''^'*''''*'- <^' 

Further integration of /i and I2 depends on the knowledge of the correlation function p^. In 
this paper p^ is assumed to have a form given by (2) . Substitute (2) in (33) and (34) and inte- 
grate with respect to x' and 7/', 

'^-''^Jy=J,, = _[-2s'(l-2772)/P.-i]"'[-22'(l-2T72)//'.-j]''^' ^ -* 
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Here the symbol Im stands for ^'the imaginary part of/^ and the wave parameters are defined by, 

D,=W{2-y')lll^. D,^4y^{z-Y)/llz, (37) 

and their significance was discussed in section 2. Since the thickness of the slab is much larger 
than the correlation distance of the irregulatities, contributions to the integrals /i and I2 come 
predominately from the region z^^lz. Therefore, (1) the limits of integration with respect to 
z^ can be extended from —00 to +00 as in (35) and (36) without introducing appreciable errors, 
and (2) because of the additional condition (17) terms like 22\l—2y^/z)/P in (35) and z''^/zp in 
(36) can be neglected as compared with unity. Introducing these approximations in (35) 
and (36), /i and I2 are respectively given by, 

h^bir^'H,, (38) 



T = r^\ M(i+z^i)^/^(i+j>i)^/^+(i-g.i9,)] , 

The exact integration of (39) is difficult to carry out since as seen in (37) D^ and Dy are functions 
of 7^ However, if the slab is thin or under conditions that I)x and Dy are almost constant and 
can be replaced by some convenient average values D^ and Dy, then, 

h=b.^'% r V(i+:p-j(i+:p-.)+(i-:p.p.) "| "' ,^^ 

Since from (26) and (27), 

<«^> = 6'<M'>(/i+/2)/2;r, (41) 

and 

<S^> = e^<M^>(/i-/2)/2T, (42) 

the mean square values are therefore given by, 

<e(0,0,.)^>=*-!<^.'%6|l+V25M(^5^Ml (43) 

<^(0,0,.)->=^^.'/-Z.6Jl-V3Sa)(l+^-.)-(i-AP.)l. (44) 

2 '1^ V 2(1+0|)(1+Z?P J 

Some special cases of interest will be considered in the following: 

(1) Z)>>1. In this case it is found that the mean square values of the phase and the 
logarithmic amplitude are equal and independent of the dimensions of the irregularities 
transverse to the direction of propagation, 

<Q'> = <S'> = e'<y?>ir"H,hl2. (45) 

This result reduces to that obtained by Chernov ^ [1961] when the irregularities have a spher- 
ically symmetric Gaussian correlation. 

(2) jD<C<C 1. In this case J9's are very small and (43) and (44) respectively become, 

<Q'> = e'<y.'>ir"H,h^ (46) 

and 

<S2>=^!<i^^ lMm+2D,Dy+Wl). (47) 



1 His eq (141) on page 75. 
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If, in addition Ix=Iv=It, hence Dj.=Dp=Drp, 

,2^ 2\ 1/2 

<gii>=^ ^ ^y^ l^bDl. (48) 

It is interesting to note that for a thin slab with spherically sj^nnietric irregularities {lx=ly= 

D^4.ac/[l\a+c)]. (49) 

Put (49) into (48) it is seen that (<^>S'^>)^^^is halved when the transmitter moves from injfinit}^ 
(a=oo) to a distance above the slab equal to that between the slab and the receiver (a=c). 
When the slab is not thin the average value of the wave parameter squared in (48) must 
be found according to the integral (39), i.e., 

=W' {I 2*[(«+*) -^*J-| 4ia+by-a*]+l [(a+by-a']y- (50) 

Equation (50) can be simplified if both the transmitter and the receiver are immersed in the 
irregularities. In this case a=0 and b = 2, substitution of (50) into (48) yields, 

<S»>=^ e'<n'>2HS'T. (51) 

It is seen that in this case the mean square values depend on the size of the irregularities 
along the direction of propagation as well as transverse to it. Equations (46) and (51) reduce 
to those obtained by Karavainikov [1957] when Z's are all equal. 

(3) D^s^l. Fortunate^ or unfortunately, in most of the ionospheric applications the 
previous approximations are invalid and hence it is necessar}^ to work with more complex 
expressions (43) and (44) of the original integral (39). Now it is seen that <C$^> and <CS^^ 
depend on D^s in a different manner. If the scales of the irregularities are determined by a 
correlation method the height of the region of irregularities can be roughly estimated by making 
simultaneous measurements of <CQ^> and <CiS'^> as was done b^^ Hewish [1952]. 



5. Transverse Correlation Functions 

The general expressions for the correlation functions were derived in section 3. In order to 
evaluate these integrals, place one receiver at (—x/2j 0, 2) and one at (x/2, 0, 2). Since the 
spacing of these receivers is perpendicular to the direction of propagation the correlations of 
the field at these two receivers are called transverse correlation functions. Approximations 
similar to those given in the last section are, 

{R,+r[)-r,=[y['+(x[-x2[/22y]/2^[, (52) 

and 

{B2+r2)-r2=[yf+(x2-x2y22y]/2^i (53) 

Define, 

I^=^Me'<^.'»{<Q,Q2>T+<S^S2>T), (54) 

and 

I,= Me'<fi'»{<QiQ2>T-<S,S2>T). (55) 
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Putting (22) and (23) into above relations these integrals are found to be, 

"~^ J J H'd'2 " L 2f ; 2E J " '''' '" 

'~ ^ }} 4f;f^ °°' L 2f; + 2E ^J '^ ^'^ "^^ 

Transform the coordinates into a relative and center of mass system, and integrate with 
respect to a' and /3' as was done in the last section, 



I 



Here the limits of integration are identical to those given b}' (33) and (34). Ju fact /^ and I^ 
reduce to /i and /2, respectively, when the spacing of the receivers x goes to zero as expected. 
For the correlation function of tlu^ form givcMi ])y (2), tlu^ int(\iiral /^ is giv(Mi aj)[)roxiniat(*ly })y, 






H!'^h-<m « 



where it is seen that the arguments of these error functions are actually tlu* magnification pa- 
rameters (discussed in 2) at the bottom and top of the shib. Integrate (57) with respect io :r/ , y\ 
and z', 

For cases in which D^s avo almost constant, (59) is approxinuitely, 

U=i.%l.^/2.) I.n{,(l + .5.)- [erf Qj;^^..? )-'^^^ Un + MY'^ )]} ' ^''^ 
Define the normalized transverse autocorrelation functions by, 

Pq{x) = <Q,Q2>tKQ'>, and Ps{x) = <S:S~2>t/<S'>, (61) 

then 

Pq(^) = (/3 + /4)/(/i+/2), (62) 

Ps{x) = {h-h)l{h-h). (63) 

For the special case in which the slab is very thin, (58) and (60) can be approximated. If, 
additionally, the Z^^s are either large or small some simple results can be obtained. These are 
considered in the following : 

(1) D's > >1. Keep only leading terms in (62) and (63) and consider the region olx^l^zja. 
The correlations of the phase and the logarithmic amplitude reduce to, 

Pq{x) = ps{x)=qx^ {—a^x^/llz^). (64) 

Hence it is seen that both correlation functions are initially Gaussian and the '^scale^^ of the 
irregular waves is a factor z/a times the *'scale^^ of the fluctuations in the refractive index 
consistent with the idea of magnification discussed in section 2. 

(2) Z)'s<^l. Again keeping only leading terms, 

pQ(x)=exp {—a?x^lllz-), (65) 
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and 



Ps{x) = 



g'x'D, 



1— 



4:a'x'D, {sD,+Dy ^, 






r~2 



ll{3Dl+3Dl+2D,D,) 



exp 






(66) 



The fluctuations in the phase is a factor z/a larger than the fluctuations in the refractive index. 
However, now the fluctuations in the ampUtude is influenced by not only the fluctuations of 
the refractive index in the direction of the spacing of these receivers (x-direction) but also, 
though weakly, transverse to it (y-direction). If additionally ly=lx, then (66) reduces to, 



Ps{x}=( 1 



2aH^ 



zHl 






(67) 



Though (66) and (67) are not Gaussian they have a correlation distance of the order of Ixz/a 
and (66) approaches Gaussian in the limit ly<^lx- 

6. Longitudinal Correlation Functions 

In this case the two receivers are placed along the direction of propagation at (o, o, z) and 
(o, o, z +A2) with a separation I\z=Z2 — Zi. Proceed in the same manner as the previous two 
sections by defining, 



and 



h={irle'<tx'»{<Q,Q,>^ + <,S,S,>,}, 



h=--{i^le'<^'»{<Q,Q.2>L-<SA>j.). 



(68) 



(69) 



It can be shown that for the spacing A^ of the order of the correlation cUstance I^ and h reduces 
to /i and I2 respectively. Therefore, the phase and the logarithmic amplitude are perfectly 
correlated in a distance A2^zl/(a-\-b), 



7. Method of Generalization 

In the foregoing sections expressions of mean square values and the autocorrelation func- 
tions are obtained by assuming that the fluctuating part of the dielectric constant has a Gaus- 
sian character given by (2). As mentioned before one of the most important features of the 
ionospheric irregularities is that they are strongly anisotropic. If p^ departs from Gaussian but 
still possesses ellipsoidal symmetr}^ some of the earlier derived expressions can be generalized. 
It is convenient to define 

P=x'lll+flll+2Vll (70) 

Introduce F{t) so that p^ (p) and F{t) are related through the Laplace transform, i.e., 

/» 00 

P„(p)=J^ Fit)expi-pt)dt. (71) 

The function F{t) can be chosen to simulate almost any correlation function with ellipsoidal 
symmetry. The eflfect of introducing (71) in the integrals I i{lx,ly,lz) , ^ = 1 to 6, is simply to 
replace them by 



r dtF{t)IJk,k,k^■ 



(72) 
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Similar generalizations can be applied to <^Q^^, <CS^y>, ps, and pq if they are evaluated 
explicitly. But now the averaging of Z>'s is not permissible and tliere does not seem to be any 
simple way out. 

By means of this device the generalization of (38), (45), and (46) is equivalent to multipli- 
cation of these respective equations by the factor, 

rmdt=^r^dp, (7..) 

Jo -yjt '\7rJo -yjp 

and the factor that goes with (48) when (50) is substituted, and with (51) is 

The generalization of the correlation functions can be done in an identical manner. If 

Ma{x)=ax/2lx, (75) 

the generalization of (64) and (65) can be shown to be 



X 



J 00 




(76) 



p^{p)p '^Hp 



The generalization of (67), when (49) is valid, is approximately 
1 

Jo --W~ '' '^' 



(77) 



In the general case it can still be proved that tlu* longitudinal correlation functions are perfect 
for Az of tlie order of 2l/{a + b). 

8. Applications to the Ionosphere 

The first clue to the existence of the irregularities in electron density was obtained by Hey^ 
Parsons, and Phillips [1946J. They discovered that the intensity of the radiation from radio 
stars fluctuate on certain occasions in the VHF band. At first it was thought that this indi- 
cated the source was variable in its power output. Subsequently spaced-receiver experiments 
by Smith [1950] and by Little and Lovell [1950] gave convincing proof that the cause of these 
fluctuations is in the ionosphere. It is now generally believed that there exist in the ionosphere 
blobs with excesses or deficiencies of electrons which scatter waves irregularly. These blobs 
are elongated along the Earth's magnetic field lines with an axial ratio of five to one [Spencer, 
1955]. Past studies of these related problems using radio stars and ground-based transmitters 
have been well summarized by Little et al. [1956], Ratcliffe [1956], and Booker [1958]. 

In the following three subsections the relation between fluctuations in refractive index 
and electron density, the scintillation indices and the autocorrelation functions are considered. 

8.1. Fluctuations in Refractive Index and Electron Density 

In case of ionospheric applications for a frequency of 20 Mc/s or higher the longitiKhnal 
approximation to the Applet on-Hartree foinuda [Appleton, 1932] is valid uidess the ray is 
within a few degrees from the direction of exactly perpendicular to the Earth's magnetic field. 
Using th(^ standai'd notations reconunended by U. R.S.I. [Ratcliffe, 1959] tlie refractive index 
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is given by, 

where AX (x) is tlie part due to electron density fluctuations and is a stochastic function of 
position. The upper sign in (78) is for the ordinary ray and the lower sign for the extraor- 
dinary ray. Taking the square root of (78) and keeping only first order terms, the refractive 
index is, 

*=<->(-2<i#xj> <™) 

Comparison of (79) with (1) shows that the expression for the mean square of the percent 
fluctuations in the refractive index as, 



where 



^<M^>=^<(AiV)^>, (80) 

=2.8178X10"^^ m (classical electron radius), 



^ 47r€omc^ 
AA^= electron density fluctuations. 

Since from (79) the fluctuations in the refractive index are proportional to the fluctuations 
in the electron density, the normalized autocorrelations of both must be identical. Therefore, 
the formulas derived in the previous sections can be applied by replacing the mean square 
electron density through the relation (80) and by noting that Vs now represent the scale of 
ionization irregularities. 

8.2. Scintillation Indices 

With the advent of artificial satellites it is now possible to carry out a series of measurements 
with a transmitter at a variable height. The transmitter may be below, within, or above 
the slab of irregularities. Hence it is desirable to know how the scintillation indices of the 
amplitude and the phase (defined as<(^^>^/^ and <^S^^'^^^ respectively) vary with the height 
of the transmitting satellite. The result is shown in figure 4. Three cases have been con- 
sidered: (1) The direction of propagation paralley to the major axis of the irregularities and; 
(2) two cases when the direction of propagation is transverse to it. The numerical values 
chosen are tabulated below: 





2XlO^-\a\ when-2X104<a<0. 




b=4 




[2X10^ whena>0, 


In all cases: 






c= 12X10^ 


Case 1 


Longitudinal case: 




1,=1,=4:X10', Z, = 2X10l 


Case 2 


Transverse case 1: 




/^-2X10^ ly=l^4:XlO'. 


Case 3 


Transverse case 2: 




Z,-L-4X10^ 1,=2X10\ 
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In terms of physical units, for A' = 400/kin (corresponding to a frequency of 19.1 Mc/s if the 
refractive index is 1 ) t hese values are given by, 

i_f^Okm—\a\ when -50 ktn <a<0, 
\50 km when a>0, 



c = 30i) km, 



Case 1 Longitudinal case: 



h=fy=^ km, h = 5 km, 



Case 2 Transverse case 1: 



L=5 km, 



ly=lz=^ km. 



Case 3 



Transverse case 2: 



l,=l,= l ki 



Ji, 



L = 5 Ivjn, 



In calculation the approxhnate fonnulas (43) and (44) have been used. These fornmlas can 
be simplified considerably in ease 1. Since Dx and Dy enter in a synnnetrical manner, the cor- 
responding scintillation indices in case 3 are identical to those in case 2. The actual plotted 
values are (2<r/>/6<M'> ^'^ '^'XUr' and {2< S'> / e< fM'> . ir'^')'^'X10-\ These nat- 
urally are proportional to the scintillation indices defined in this paper. As an a])|)roximation 
Z>'s are calculated from (37) by assuming 7' at the midway of the slab. When the satelhte is 
inside the slab the effective slab thickness is reduced accordingly. 




CASES 283 



10x10 (NORMALIZED) 



-50 100 200 

HEIGHT OF TRANSMITTER ABOVE TOP OF SLAB, a 



Figure 4. Dependence of scintillation on the height 
of the transmitter. 

By careful examination of figure 4 the following interesting points are noted : 

(1) The scintillation index of the amplitude varies monotonically as a function of the 

height while the scintillation index of the phase rises rapidly when the transmitter is within 

the slab and reaches a maximum value near the top of the slab and then decays slightly above 

the slab. Therefore, it appears that the measurement of phase scintillation may enable us 
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to determine the height of the bottom, the thickness, and the height of the top of the slab if it 
exists at all. The existing experimental data are rather meager. Some rough estimates 
using signals from satellite 1958 A2 on 20 Mc/s have been made [Swenson and Yeh, 1961]. 
The amplitude scintillation roughly agrees with figure 4. (See their figure 6, which is reproduced 
here as figure 5. Note that the scintillation index defuied in that paper is identical to that 
defined here only in the limit of very weak amplitude fluctuations.) The data presented in 
figure 5 were collected in Baker Lake, Canada which is north of the auroral zone. There is 
evidence that the region of irregularities may var}^ in height by more than 500 km [Easier 
and DeWitt, 1962] in the auroral zone. This may account for the spread of these points. 
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Figure 5. Experimental results of scintillation de- 
pendence on the height of the satellite. 

(2) For a constant thickness of the slab figure 4 shows that there is more scintillation in 
the longitudinal case (case 1) than that in the transverse case (case 2 and case 3). The phase 
scintillation is 2 to 2.2 times more and the amplitude scintillation is 3.1 to 3.7 times more in 
the longitudinal case than the corresponding values in the transverse case. In the actual 
ionospheric observations the effective thickness varies according to the secant of the zenith 
angle of the 'direction of propagation. Since these irregularities are alined along the magnetic 
field lines and the dip angle at, for example, the University of Illinois is 71.3°, the effective 
thickness in the longitudinal case (case 1) is then a factor tan 71.3° = 2,95 smaller than that 
in the transverse case (cases 2 and 3). This will bring curves of cases 2 and 3 closer to the 
corresponding curves of case 1. Hence in actual ionospheric observations at temperate and 
high latitudes the scintillations should not change appreciably with the zenith angle of the 
transmitter and preliminary results of the observation of satellite signals at various stations 
operated by the University of Illinois seem to indicate that this is so. However, for a station 
near the magnetic equator we would expect that the effect of the effective thickness nia}^ 
enhance the differences in scintillations considered above and the signal from an orbiting 
satellite at low angles to the north or south of the station should exhibit stronger scintillations 
than when it is overhead. 

8.3. Transverse Correlations 

It has been shown in section 6 that the longitudinal correlations are perfect for a separa- 
tion A^ of the order ^//(a + 6). Therefore only the transverse correlation will be considered 
here. 

Examining the integral I^ given by (59) or its approximate expression (60), it is seen that 
the correlation functions depend mainly on the dimension of the irregularities along the spacing 
of the two receivers and not on their dimension transverse to the spacing only if the wave 
parameter is small compared to 1 . That this is intuitively true can perhaps be argued physi- 
cally as follows. When Z)<^<^1 the irregularities are large compared with the first Fresnel 
zone. Since most of the contributions to the signal at the receiver come from the scattering 
in the first Fresnel zone it seems reasonable to expect that the correlations should not be sensi- 
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tive to the dimension of the irregularitic^s transverse to tlie spacing of the receivers. However, 
when i9>> 1 tlie irregnhirit ies are small compared with the first Fresnel zone and the scattering 
is all important toward (he end of (he irregularity. Therefore, in this case the transverse 
dimension of the irregularides a(r(M'(s the correlation function. 

As specific nmneiMcal e\ajnj)l(^s the three cases discussed in the last subsection are again 
considered. The numerical values are identical with those tabulated previously except that 
now a fixed value of a= 12X10* is taken. For ^=400/km this is equivalent to a height of 300 
km above the top of the slab. The curves are shown in figure 6. The ordinate is tlu^ inagni- 
fication parameter defined in (7) where 7' is assumed to be the distance from tin* transniitler 
to the center of the slab, i.e.. 



M,-- 



{a+b/2)x 



(81) 



In calculating these curves formulas (58) and (60) have been used. 

Experimentally one usually obtains the fading records at a fixed station froju a moving 
transmitting vehicle. For all practical purposes it is equivalent to the case considered in this 
paper. The autocorrelation function of the amplitude obtained experimentally (defined in a 
slightly different manner) agree reasonably well with those theoretical curves shown in figure 
6 [Kent, 1959; Swenson and Yeli, 1961]. 
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a = 12 X 10 

b = 2 X 10 
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Figure 6. Transverse correlation functions. 
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